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ABSTRACT 

The cusped NFW universal density profile suggested by typical CDM models has 
been challenged in recent years by the discoveries of the soft cores with finite central 
density for a broad range of masses from dwarf galaxies to clusters of galaxies. It is thus 
desirable that a new, analytic model would instead become available for virialized dark 
halos. One promising candidate is probably the empirical density profile proposed by 
Burkert (1995), which resembles an isothermal profile with a constant core in the inner 
region and matches the NFW profile at large radii. Meanwhile, such a revised dark 
halo (RDH) profile has turned out a great success on galactic scales. This stimulates 
us to apply the RDH profile to more massive systems like clusters of galaxies. In this 
paper we have made an attempt to derive the radial density profile of intracluster gas 
from the RDH profile under the isothermal and hydrostatic equilibrium hypotheses, 
and compare it with those revealed by X-ray observations and inferred from the NFW 
profile. It is shown that the RDH predicted gas density can be well represented by 
the conventional (3 model with a typical (3 parameter of (3 ~ 0.7-0.9. Alternatively, 
fitting the theoretically predicted X-ray surface brightness profile to an ensemble of 45 
X-ray clusters observed by ROSAT, we find that the RDH and NFW profiles become 
to be almost indistinguishable from each other, and their characteristic density and 
scale length parameters are strongly correlated. Yet, unlike the NFW model, the RDH 
profile can allow us to work out straightforwardly the central dark matter density 
from X-ray measurements of the surface brightness and temperature of clusters. It 
appears that the resulted central densities of the 45 clusters have an average value of 
(po) ~ 0-01 M pc~ 3 , in agreement with the result estimated on galactic scales, which 
reinforces the claim for the presence of the soft halo cores over the entire mass range. 

Subject headings: cosmology: theory — dark matter — galaxies: clusters: general - 
X-rays: galaxies 
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1. Introduction 

For decades many efforts have been made towards the understanding of the radiai density 
profile of intracluster gas from the well-motivated physical mechanism, in an attempt to recover 
the empirical (3 model (Cavaliere k Fusco-Femiano 1976) which fits nicely the X-ray observed 
surface brightness distribution of clusters. Assuming both galaxies and gas are the tracers of the 
shape and depth of a common gravitational potential of a cluster, Cavaliere k Fusco-Femiano 
(1976, 1978) obtained an analytic gas density profile resembling the (5 model in shape if the King 
model is used to represent the galaxy number density profile in the inner cluster region. In recent 
years, the rapid progress of numerical simulating techniques has permitted the reconstruction of 
the dark matter halos with an unprecedented resolution, ranging from galactic scales of ~ 1 kpc to 
large-scale structures of ~ 10 Mpc. This leads one to view the issue at a different angle: Given the 
gravitational potential wells defined by the dark halos, how is the intracluster gas distributed in 
clusters if the hydrostatic equilibrium between the gas and the underlying gravitational potentials 
has been built up ? In particular, much attention has been paid to the issue of how to reconstruct 
the radial gas density and temperature profiles if the dark halos follow the so-called universal 
density profile (Navarro, Frenk k White 1995; NFW) suggested by high-resolution simulations 
within the framework of typical CDM models such as SCDM, LCDM and ACDM (Makino, Sasaki 
k Suto 1998; Suto, Sasaki k Makino 1998; Yoshikawa k Suto 1999; Wu k Chiueh 2000). It 
turns out that there is indeed a striking similarity between the predicted X-ray surface brightness 
profiles of clusters and the conventional (5 model. This has stimulated several authors to apply the 
NFW predicted X-ray surface brightness profiles to the observed ones for an ensemble of X-ray 
clusters (Makino k Asano 1999; Ettori k Fabian 1999; Wu k Xue 2000; Wu 2000). 

Yet, besides its uncomfortable singularity at r = 0, the cusped NFW profile has been shown 
to be in conflict with observations (e.g. Tyson, Kochanski k Dell'Antonio 1998; Navarro k 
Steinmetz 2000; Firmani et al. 2000; and references therein). In fact, it has been noticed that, 
while the NFW profile yields a gas density profile close to the empirical (3 model, the core radius 
predicted by the cusped density profile may be smaller than the actually observed one (Makino 
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et al. 1998). Moreover, the NFW profile leads to an increasing gas temperature towards cluster 
centers (Wu & Chiueh 2000), in contrast with the presence of the cold gas components detected 
very often inside the X-ray cores. Motivated by the soft inner matter distributions of the dark 
halos revealed observationally from dwarf galaxies to rich clusters, Spergel & Steinhardt (2000) 
recently proposed that the CDM particles are self-interacting. As a result, the collisional CDM 
particles, in a similar way to the baryonic particles, will produce the less centrally concentrated 
structures. This weakly self-interacting CDM model has soon attracted many investigations, 
among which the numerical simulations of structure formation based on some simple physical 
consideration of the self-interacting CDM particles have successfully provided a scenario that is 
essentially consistent with the existing observations (Hannestad 1999; Burkert 2000 and references 
therein). 

As a natural extension, one may address the following question: Does there exist a similar 
analytic expression to the NFW universal density profile that one can use to approximately describe 
the matter distribution of the virialized dark halos resulted from the weakly self-interacting CDM 
model ? A conclusive answer to such a question seems not easy: Unlike the standard CDM 
particles, for which there is no need to consider the interaction between particles except their 
gravity, the collisional, warm CDM particles contain an unknown parameter, the cross-section. 
Consequently, even with the help of high-resolution simulations, it is still hard to completely 
determine the final configurations of dark halos. Under present circumstances, one promising 
candidate is probably the empirical density profile suggested by Burkert (1995): 

pDMir) = - — ■ 2 (1 

where po and tq are the central density and the scale length, respectively. This revised dark 
halo (RDH) density law resembles an isothermal profile in the inner region with a constant core 
vq, while in the outer region the mass profile diverges logarithmically with r, in agreement with 
the NFW profile. So, such a density profile does have the desired properties at the central and 
outermost regions. In particular, it fits fairly well the dark matter distributions of dwarf galaxies 
revealed by both the rotation curves and the numerical simulations of evolution of halos consisting 
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of weakly self- interacting CDM particles (Burkert 2000; Salucci & Burkert 2000). 

Motivated by the apparent success of the RDH profile of eq.(l) on galactic scales and the 
possible existence of weak interaction between CDM particles, in this paper we would like to apply 
the RDH profile to more massive systems like clusters of galaxies. We would like to demonstrate 
the radial density profile of intracluster gas tracing the gravitational potential defined by the RDH 
profile, and compare the expected X-ray surface brightness with the X-ray observations and other 
models (e.g. the conventional (3 model and that predicted by the NFW profile), although we have 
an intuition that the RDH profile, as a combination of the NFW profile and the (3 model, would 
provide an essentially similar result. This study will nevertheless constitute an important test for 
the universality of the RDH profile as the virialized dark halos over the entire mass range, and will 
also allow us to examine whether there are any common properties in the NFW and RDH profiles. 
Eventually, it is hoped that this study will be useful for a conclusive answer to the question as to 
whether the RDH profile can be used to replace the role of the NFW profile for the virialized dark 
halos, if CDM particles are indeed weakly self-interacting. 

2. Density profile of intracluster gas 

2.1. Self-gravity of the gas: excluded 

If the dark halo of a cluster follows the RDH profile described by eq.(l), the total dark matter 
enclosed within radius r is 



where x = r/rQ. Assuming that the intracluster gas is isothermal (with a temperature of T) and in 
hydrostatic equilibrium with the underlying gravitational potential dominated by M^m, we have 



M DM {x) = 4irp rorh(x); 
m(x) = - ln(l + x) + - In (1 + x 2 ) — arctan(x) , 



(2) 



(3) 



GM DM (x) 
x 2 



jj,m p n gas (x) dx 



kTr dn gas (x) 



(4) 
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in which n gas (x) is the gas number density and \x is the average molecular weight. Note that we 
have neglected the self- gravity of the gas for the moment. A straightforward computation yields 
an analytic form of gas density profile: 



n g as(0) 

where 



e -( 1 +^)^^(l + x)( 1+ ^(l + x 2 )5^- 1 )l 2 , (5) 



4irG[im p p r% 

«o = ^ • (6) 

In order to avoid the divergence of the resulting X-ray surface brightness to be discussed below 
because of n gas (oo) = n 9 as(0)e _7ra °/ 4 , in a similar way to the treatment of the NFW predicted 
gas density profile (see Wu 2000), we introduce a normalized, background subtracted gas number 
density h gas (x) = [n gas (x) - n gas (oo)]/[n gas (0) - n gas (oo)], which reads 

1 



fbgas (#) 



7TQIQ 

— 1 



ei -(l+I)arctai.x (1 + + ^1(1-1)] ? _ j j (?) 

Another way to deal with the non-zero background gas density predicted by the dark halo model 
is to truncate the cluster at a certain radius (e.g. the virial radius) as adopted by Makino et al. 
(1998). It is indeed unfortunate that one has to add an arbitrary, unphysical constraint on the gas 
density profile to ensure the convergence of the X-ray surface brightness. 

In Fig.l we demonstrate the radial profiles of the scaled gas density n gas (x) for typical 
clusters with ao = 5, 10 and 20. A glimpse of Fig.l seems to suggest that all the resulted 
profiles of h gas {x) resemble the conventional (3 models in shape. We then overlap the (3 model, 
n* as (x) = [1 + {x/xc) 2 }-^/ 2 , to each curve with (/3, r c /r ) = (0.40, 0.86), (0.74, 0.85) and (1.56, 
0.92) for ao =5, 10 and 20, respectively, where x c = r c /ro is the scaled core radius. Yet, like 
the actual fitting of the (3 model to the X-ray observed surface brightness profile, the best-fit (3 
parameters depend also on the extension of the fitting regions. For a typical cluster of ao ~ 10, the 
best-fit (3 value over a region out to 10r c is (3 ~ 0.7, in good agreement with X-ray observations. 
In particular, the gas core radius takes roughly the value of tq. Also plotted in Fig.l are the 
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corresponding gas densities predicted by the NFW profile with a = (Makino et al. 1998; Wu 
2000) 

, . (l+x) a / x - 1 

h 9as(x) = V J_ l , (8) 

where x = r/r s and a = AmG\im p p s r^jkT . A visual examination of Fig.l reveals that there are 
some differences between the expected radial variation of intracluster gas from the RDH profile 
and that from the NFW profile, which are reflected not only by the significantly different scale 
lengths but also by the different shape. Indeed, it is unlikely that the two types of density profiles 
can be made to be identical simply by a horizontal replacement. Whether these differences are 
significantly important will be discussed when the two density profiles are both applicable to an 
ensemble of X-ray clusters (section 3) . 

EDITOR: PLACE FIGURE | HERE. 



2.2. Self-gravity of the gas: included 

The total mass in gas within radius x is simply 

Mgas(x) = AiT^m p rl \ n gas (x)x 2 dx. (9) 



o 



When the self-gravity of the gas is included, the hydrostatic equation becomes 

G[M DM (x) + M gas (x)\ _ kTr dn gas (x) 



x 



2 



Hm p n gas (x) dx 



dfb _ (1 ~ fb)\K l ~ fb)n gas - fbPD M]x 
dx m 



2 



(10) 



If we introduce the volume-averaged (gas) baryon fraction fb(x) as a new variable: 

fh(x) - 

Jb[x) M DM (x) + M gas ( x y [Ll) 

we can obtain the following two first-order differential equations 

dflqas "ff^qas /in\ 

dx (1 - } b )x z 



(13) 
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where n gas = n gas (x)/n gas (0), b = pm p n gas (0)/ p and p DM = Pdm/po = (1 + x) 1 (1 + x 2 ) 1 . We 
need to specify the boundary conditions in order to solve the above equations. The first condition 
is obviously 

n sas (0) = f. (14) 

Following the recent work of Wu & Chiueh (2000), we choose the boundary condition of fb(x) 
such that the baryon fraction within the virial radius r v i r (or c = r v i r /ro) should asymptotically 
approach the universal value fb,BBN = Qb/^M, where fi& and Qm are, respectively, the baryon 
and total mass densities of the Universe in units of the critical density p c for closure, namely, 

fb(c) = fb,BBNi (15) 

fl„ c = 0. (16) 
Here the scaled virial radius or the so-called concentration parameter c is defined by 

M DM (c) = —r% Pc c 3 A c , (17) 



(18) 



or 

m(c) A c 1 

= TV 

in which S c = po/p c , and A c represents the overdensity parameter of dark matter with respect to 
the average background value p c and will be taken to be 200 in our computation below. 

The free parameters involved in eqs.(12) and (13) and the boundary conditions are ceo, b, 
fb,BBN and c or 5 C . However, there are only two independent parameters with the restrictions of 
eqs.(15) and (16), which we choose to be 5 C and fb b b n below. We will perform the numerical 
searches for the solutions of eqs.(12) and (13) under the boundary conditions of eqs.(14)-(16). 
Technically, we search for the solutions over a two-parameter space (uq, b) by iterations until the 
boundary conditions are satisfied, which enables us to work out the radial profiles of gas density 
and baryon fraction, together with a unique determination of the parameters ao and b. In Fig.2 we 
demonstrate a set of solutions for two typical choices of fb,BBN and 5 C : fb,BBN = (0.05, 0.1) and 
5 C = (10 4 , 10 5 ), and the resulting ao and b values are listed in Table 1. It appears that the derived 
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density profiles of the baryon fraction exhibit no dramatic variation over the whole clusters: 
For a small value of S c ~ 10 4 , /& increases slightly with outward radius and eventually matches 
the background value at virial radius, while for a large 5 C ~ 10 5 , fb would reach a maximum 
before approaching asymptotically the universal value. The parameter ao also acts roughly like 
a constant (~ 12) even if the characteristic density 5 C changes by a decade, and the parameter 
b turns out to be an increasing function of fb,BBN, which arises simply from b oc fim p n gas (0). 
Now, we concentrate on the derived profiles of gas density. All the predicted density curves of 
intracluster gas shown in Fig. 2 seem to well resemble the (3 models, which is illustrated by our 
superimposed (3 model onto each of the derived density profiles. The most remarkable feature 
is that the corresponding values of (3 and core radius have all fallen into very narrow ranges of 
0.87 < [3 < 0.98 and 0.65 < r c /ro < 0.88 for our choices of the universal baryon fraction fb,BBN 
and the characteristic density 5 C for typical clusters (see Table 1). An increase of 6 C up to 10 6 , 
the roughly largest density for clusters (see next section), only leads to a minor modification to 
these limits. This indicates that the isothermal gas in different clusters should essentially follow 
a similar distribution, i.e., the observationally determined (3 among different clusters should not 
show a large scatter around the mean value (3 ~ 0.9. Note that the gas core radius r c may vary 
substantially because of the different tq for different clusters. While the predicted density profile 
of intracluster gas demonstrates a property basically consistent with what has been known for 
X-ray clusters, the theoretically expected (3 parameter is likely to slightly exceed the presently 
determined one from X-ray observation, (3 ~ 0.7. The former arises mainly from our restriction 
that the baryon fraction should asymptotically match the background value. Previous studies with 
this constraint have also arrived at a similar conclusion: the (3 value in the (3 model for intracluster 
gas is required to be larger than a certain low-limit (Wu & Chiueh 2000). Actually, it is not 
impossible that the presently fitted (3 parameters based on the X-ray observed surface brightness 
profiles of clusters are biased low because of the influence of the cooling-flows and the small fitting 
regions. Excluding the cooling flows or adopting a double (3 model fit may moderately raise the 
observationally determined (3 values, giving rise to (3 0.7-0.8 (Vikhlinin et al. 1999; Xue & Wu 
2000). Of course, the present computation has been made within the framework of isothermality, 
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and the above prediction cannot be taken too literally unless the non-isothermal gas is included. 

EDITOR: PLACE FIGURE | HERE. 

Considering the fact that fb exhibits only a minor variation over the whole cluster, and also 
for illustrating the effect of the self-gravity of the gas, we can provide an approximate and analytic 
form of the gas density by taking fb = fb,BBN in eqs.(12) and (13). Consequently, eq.(12) reduces 
to eq.(4) and the gas density is given by the analytic expression eq.(7), in which the parameter 
«o is now replaced by ag = ao/(l — fb)- In Table 1, we list the f3 parameters and core radii in 
the (5 model fits to our approximate solutions using the same input values of fb,BBN and 5 C (or 
ao). In the case of the low density 5 C = 10 , the agreement between the exact and approximate 
solutions is fairly good, while the approximate solutions seem to underestimate the (5 values by 
~ 0.15 for S c = 10 5 . The latter is partially due to the fact that the fitting of the f3 model is 
made over rather a large region out to x = c = 14.2, where a shallower density profile occurs 
according to the RDH profile (see Fig.l). Another reason is that in the case of 5 C = 10 4 , the 
baryon fraction remains roughly unchanged across the clusters with a relative variation rate of 
\fb(c) - f b (0)\/f b ,BBN < 20%, in comparison with |/ 6 (c) - f b (0)\/ f b , BBN ~ 100% for 5 C = 10 5 . 
Namely, the small 5 C clusters seem to meet more easily the condition for the approximate solutions 
(/ft = fb,BBN) than the large 5 C clusters do. Nevertheless, we conclude that the derived gas 
distributions with and without the inclusion of the self-gravity of the gas do not show very 
significant difference, and in general, the self-gravity of the gas only leads to a slightly steeper gas 
density profile as a result of the increase of the underlying gravitational mass. This is consistent 
with the similar study for the NFW profile (Suto et al. 1998). 

3. Application to X-ray clusters 

In this section we conduct a comparison between our derived gas density profile and X-ray 
observations, which are linked up through the X-ray surface brightness profiles of clusters, S x . In 
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the scenario of the optically thin, isothermal plasma emission, 

roo 

S x (x) cx / h 2 gas d£, (19) 

Jx 

where the integral is performed along the line of sight £. Here we take the analytic form of the gas 
density eq.(7) and neglect the contribution of the gas self-gravity. We intend to fit our theoretically 
expected X-ray surface brightness profile eq.(19) to an ensemble of the X-ray observed surface 
brightness profiles of clusters, which will allow us to determine the characteristic parameters, po 
and ro, for the RDH profile and compare them with the results from other models, e.g. the (3 
model and the NFW profile. 

We use the ROSAT PSPC observed surface brightness profiles of 45 nearby clusters compiled 
by Mohr, Mathiesen & Evrard (1999, MME). Several models have been already tested with this 
sample, such as the (3 model, the double (3 model and the NFW predicted density profile (MME; 
Wu & Xue 2000; Xue & Wu 2000). In a similar way to our previous analysis for the NFW profile 
(Wu & Xue 2000), we perform the x 2 -fit to get the best-fit parameters ao and ro, in which we keep 
the same outer radii of the fitting regions as those defined by MME. In order to examine how our 
results are affected by the presence of the cooling flows in some clusters, we perform our fittings by 
using the entire data points of S x (model A) and excising the central region of 0.05 Mpc in each 
cluster (model B), respectively. For the latter the reason that we adopt the same inner radius for 
all the clusters is to guarantee the uniformity of the excision (Markevitch 1998). We then compute 
po in terms of eq.(6) by taking the X-ray temperature data from the literature (see Wu, Xue & 
Fang 1999; and references therein). For majority of the clusters we use the cooling flow corrected 
temperature data by White (2000). In Fig.3 we illustrate a typical example of the observed and 
our fitted surface brightness profiles (model A) for cluster A3158. For comparison, we have also 
plotted the results of the (3 model and the NFW profile. Essentially, these three models provide 
more or less an equal goodness of fit to the observed data with xt = 1-12, 1.24 and 1.10 for the (3 
model, the NFW and RDH profiles, respectively. From the fitting of the X-ray surface brightness 
profile alone, it may be hard to reject any of these models. In Table 2 we list the best-fit values 
of po and ro for the 45 MME clusters by model A, together with the results for the (3 model and 
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the NFW profile, in which all the quoted errors are 68% confidence limits. The Hubble constant 
is taken to be Hq = 50 km s" 1 Mpc" 1 . When the central regions of 0.05 Mpc are excised in the 
fittings, the best-fit parameters for all the three models are only moderately affected. However, 
this can significantly improve the goodness of the fittings characterized by the reduced xt- After 
the excision of the central regions (model B), the fractions of clusters with <2 increase from 
(10/45, 14/45, 7/45) to (26/45, 25/45, 22/45) for (/3, NFW, RDH) models, respectively. Recall 
that the similar fraction of 26/45 is found by MME for the single (3 and double f3 model fittings. 
It appears that about half of the clusters cannot be well fitted by any of these models even if the 
central regions are excised. 

EDITOR: PLACE FIGURE | HERE. 



3.1. RDH vs. (3 model 

We first analyze the possible link between the RDH profile and the (3 model. For this purpose 
we plot in Fig. 4 the best-fit ao vs - P an d tq vs. r c for the 45 MME clusters obtained by model A. 
It is immediate that there exist strong correlations between the scale and slope parameters in the 
two models. We fit these correlations to a power-law function, which reads 

ao = 1 1-08±0-01/?0-67±0.06 ; (A) . (2Q) 

a = 10 i.io±o.oi /3 o.86±o.04 j (B)> (21) 

and 

r = i -oo7i±o.oi5 r o.75±o.o 2) (A) . (22) 

ro = l -°-°33±0.012 r 0.83±0.Cn (B) _ (23) 

Meanwhile, applying a linear fit to the data set yields 

a = (14.11 ± 1.09)13, (A); (24) 

a = (13.44 ± 0.55)0, (B), (25) 
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and 

ro = (1.46 ± 0.41)r c , (A); (26) 
r = (1.30 ± 0.18)r c , (B). (27) 

It appears that the resultant relationships with and without the excision of the central regions in 
the fits of the X-ray observed surface brightness profiles are roughly consistent with each other. 

EDITOR: PLACE FIGURE ED HERE. 



3.2. NFW vs. p model 

For comparison we display in Fig. 5 the correlations between the corresponding parameters in 
the NFW profile and the (5 model. The best-fit a-f3 and r s -r c relations using all the 45 data points 



arc 



a = 10 1.28±0.02 /3 1.46±0.08 ) (A) . (2g) 
a = 10 1.36±0.02 /3 1.81±0.08 ) ^ (2g) 



and 



r s = io«-^±u-^ r r.u a± u. Ua) (A) . (3Q) 



->0.77±0.02„1.27±0.02 



(B). (31) 



However, there are four clusters, A119, A1367, A1656 and A2255 showing a large dispersion on 
the a-(3 plane, which may significantly affect the above fittings. The best-fit a-j3 relation without 
these four clusters becomes 

a = 1O 1.16±O.O101.O1±O.(B (A); (32) 
a = 10 1.41±0.06 /3 1.26±0.01 ) (B) _ (33) 
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Alternatively, the linear fit to the data set of 45 clusters gives 

a = (15.40 ± 2.83)/?, (A); (34) 

a = (15.23 ±2.54)/?, (B), (35) 

and 

r s = (4.17±1.19)r c , (A); (36) 

r s = (3.70 ± 1.06)r C! (B). (37) 

Within the uncertainties, the last four relations are consistent with the findings by Ettori & Fabian 
(1999) based on 36 high-luminosity clusters: a = 14.34/? and r s = 3.17r c . 

EDITOR: PLACE FIGURE | HERE. 



3.3. RDH vs. NFW 

We now compare the RDH and NFW profiles. Fig. 6 shows the correlations between the 
density and length scale parameters in the two models. Applying the x 2_ fit with the inclusion of 
the measurement uncertainties to the data set gives 

_ \ / „ \ 0.73±0.02 

P0\ _ 10 1.97±0.06 / M ; (A); (38) 



and 



Pel \Pc 

Pi 

Pc 



/ n \ 0.66±0.01 

P°\ ^ 10 2.18±0.06 ^\ ; (B)) (39) 



ro = l -0-^±0.01 r 0.71±0.02 ; (A) . (4Q) 

ro = l -O.S4±0.01 r 0.64±0.03 ) (B) _ (41) 

The average ratios of po/p s and ro/r s are, respectively, 

po = (8.63 ± 5.61)p s , (A); (42) 

p = (7.89±4.93)p s , (B), (43) 
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and 



r = (0.36 ± 0.09)r. 



(A); 



(44) 



r = (0.37 ± 0.10)r. 



(B). 



(45) 



The last two relations are consistent with the estimate of Burkert (2000), r$ ~ 0.2r s , for 
dwarf galaxies. The existence of these strongly positive correlations is helpful for us to make a 
quantitative comparison between the two models. Indeed, from their predicted gas densities, and 
in turn their X-ray surface brightness profiles of clusters, along with the correlations between the 
characteristic density and scale length parameters, we are unable to distinguish the two models as 
the dark halos of clusters. However, there is a remarkable difference, the central density. In fact, 
the characteristic density parameters, po an d p s , in these two models have very different meanings. 
po represents explicitly the central density of dark matter, while the parameter p s corresponds to 
the density nowhere in clusters. Note that the large error bars in the linear po~p s relation (eqs.[42] 
and [43]) and the average values of p s and po listed in Table 2 are mainly due to the inclusion of a 
very few clusters (e.g. A262 and A3526) whose X-ray surface brightness profiles show two distinct 
length scales 



Recall that the disagreement between the shallower central density profiles required by 
various observations with the cusped density profile provided by the NFW model is one of the 
primary motivations for advocating the scenario of the weakly interacting dark matter particles 
(Spergel & Steinhardt 2000). A combination of the RDH profile and the X-ray surface brightness 
measurements of clusters can now allow us to determine straightforwardly the central dark matter 
densities (po) of clusters. It is easy to show that pq is related to the central total mass density pp_ 
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3.4. Central dark matter density 
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inferred from the conventional isothermal /3 model through 

^ = ^G^l = 9 UJ [rj P °- (46) 
Using the linear relations established above between a > Pi r o an d r c> we have pp$ ~ po- Of course, 
the good agreement between pp a and po could be interpreted as the consequence of the common 
working hypothesis behind the two models: The intracluster gas is assumed to be isothermal and 
in hydrostatic equilibrium with the underlying gravitational potential. In a recent work, Firmani 
et al. (2000) have demonstrated an essentially constant central density of dark halos for a broad 
range of masses from dwarf galaxies to clusters of galaxies. However, there are only three data 
points on cluster scales which are derived from gravitational lensing. Now, we superimpose our 
derived central densities of dark halos of the 45 MME clusters on their illustration of halo central 
density vs. maximum rotation velocity (Fig. 7), in which the velocity dispersion is plotted as the 
horizontal axis for clusters. Although the central densities of the clusters span almost two decades 
from 10 -3 to 10 _1 M pc~ 3 , the average value (p ) = 0.012 M Q pc~ 3 (O.OO6M pc~ 3 for model 
B) is in agreement with the one (0.02 M Q pc -3 ) reported by Firmani et al. (2000). This reinforces 
the claim for the presence of the soft halo cores over the entire mass range. 

EDITOR: PLACE FIGURE g HERE. 



4. Discussion and conclusions 

None of the present numerical simulations based on typical CDM models can reproduce the 
soft cores of dark halos detected observationally in various systems from dwarf and low surface 
brightness galaxies to clusters of galaxies. This challenges the analytic and elegant form of the 
universal density profile suggested by NFW as the virialized dark halos. Without a sophisticated 
treatment of the dynamical evolution of the CDM particles, we have made an attempt to adopt 
the empirical density profile (RDH) proposed by Burkert (1995) to replace the cusped NFW 
profile for dark halos. The RDH profile with a finite core rg has been shown to give a perfect fit 
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to the observed rotation curves of dwarf galaxies. In particular, it also provides an excellent fit to 
the dark matter distributions revealed by numerical simulations (Burkert 2000) in the scenario 
that the CDM particles are weakly self-interacting, suggested recently to overcome the difficulties 
of CDM models (Spergel & Steinhardt 2000). In the present paper we have applied the RDH 
profile to clusters of galaxies and derived the radial distribution of intracluster gas under the 
isothermal and hydrostatic equilibrium hypotheses. It turns out that the RDH resulted gas density 
resembles the conventional (3 model, although a slightly large (3 parameter ((3 ~ 0.7-0.9) may be 
required. The latter nevertheless agrees with the results usually found by numerical simulations 
(e.g. NFW; Eke, Navarro & Frenk 1998), and is also marginally consistent with the (3 model 
fit by excising the cooling flow regions (e.g. Vikhlinin et al. 1999) or adopting an additional (3 
model for the central X-ray emission (Xue & Wu 2000). Except for the very differently asymptotic 
behaviors at r — > 0, the RDH and NFW profiles become to be almost indistinguishable from 
each other within the extent of the current X-ray observations of clusters. By fitting the RDH 
predicted X-ray surface brightness profile to the observed one for an ensemble of 45 clusters, we 
have estimated the typical central density po and core radius tq in the RDH profile, which show 
a fairly strong correlation with the characteristic density p s and scale parameter r s in the NFW 
profile, respectively. Meanwhile, our derived central dark matter densities of the 45 clusters have 
an average value of (po) ~ 0.01 M Q pc -3 , in agreement with the result estimated on galactic scales 
(Firmani et al. 2000). 

Yes, a conclusive justification for whether the RDH profile can be used as a good approximation 
of dark halos and therefore, replace the role of the NFW profile for typical CDM models will 
be provided by high-resolution simulations incorporated with new physical mechanism for dark 
matter particles such as the weak self-interaction. Several new models with the help of numerical 
simulations are being constructed by a number of authors. It can be predicted that a new analytic 
model with an isothermal core, if it is not the RDH profile, for the structure of virialized dark 
halos will soon be available. The rotation curves of dwarf galaxies, the total mass distribution of 
clusters revealed by gravitational lensing, and the X-ray properties of clusters as explored in the 
present paper will constitute a critical test for the new RDH profile. 
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Table 1. Numerical results for typical clusters 



fb,BBN 


5 C 


c 


a 


b 


P 


r c /r 


13* 


(r c /ro)* 


0.05 


10 4 


5.48 


12.0 


0.05 


0.867 


0.705 


0.90 


0.66 


0.05 


10 5 


14.2 


13.7 


0.09 


0.974 


0.681 


0.83 


0.68 


0.10 


10 4 


5.48 


10.8 


0.09 


0.927 


0.884 


0.90 


0.63 


0.10 


10 5 


14.2 


12.8 


0.23 


0.981 


0.649 


0.83 


0.68 



*The (3 model fit to the approximate solution eq.(7) by assuming = 
fb,BBN and replacing a by a /(l - fb,BBN)- 



20 



REFERENCES 



Burkert, A., 1995, ApJ, 447, L25 

Burkert, A., 2000, ApJ, submitted (|astro-ph/0002409D 



Cavaliere, A., & Fusco-Femiano, R., 1976, A&A, 49, 137 
Cavaliere, A., & Fusco-Femiano, R., 1978, A&A, 70, 677 
Eke, V. R., Navarro, J. F., & Frenk, C. S., 1998, ApJ, 503, 569 
Ettori, S., & Fabian, A. C, 1999, MNRAS, 305, 834 

Firmani, C, D'Onghia, E., Avila-Reese, V., Chincarini, G., & Hernandez, X., 2000, MNRAS, in 
press 

Hannestad, S., 1999, preprint astro-ph/991255£ 



Makino, N., & Asano, K., 1999, ApJ, 512, 9 

Makino, N., Sasaki, S., & Suto, Y., 1998, ApJ, 497, 555 

Markevitch, M., 1998, ApJ, 503, 77 

Mohr, J. J., Mathiesen, B., & Evrard, A. E., 1999, ApJ, 517, 627 (MME) 

Navarro, J. F., Frenk, C. S., & White, S. D. M., 1995, MNRAS, 275, 720 (NFW) 

Navarro, J., F., & Steinmetz, M., 2000, ApJ, 528, 607 

Salucci, P., & Burkert, A., 2000, ApJ, in press 

Spergel, D. N., & Steinhardt, P. J., 2000, Phys. Rev. Lett., in press 

Suto, Y., Sasaki, S., & Makino, N., 1998, ApJ, 509, 544 

Tyson, J. A., Kochanski, G. P., & DelPAntonio, I. P., 1998, ApJ, 498, L107 

Vikhlinin, A., Forman, W., & Jones, C, 1999, ApJ, 525, 47 

White, D. A., 2000, MNRAS, 312, 663 

Wu, X.-P, 2000, MNRAS, in press 



- 21 - 

Wu, X.-R, & Chiueh, T., 2000, ApJ, submitted 

Wu, X.-R, k Xue, Y.-J., 2000, ApJ, 529, L5 

Wu, X.-R, Xue, Y.-J., & Fang, L.-Z. 1999, ApJ, 524, 22 

Xue, Y.-J., & Wu, X.-R, 2000, MNRAS, submitted 

Yoshikawa, K., & Suto Y., 1999, ApJ, 513, 549 



This manuscript was prepared with the AAS MpX macros v4.0. 



-22- 



Fig. 1. — The radial density profiles of intracluster gas derived from the RDH profiles (solid lines) 
and the superimposed ft models (open squares). The parameters (a^, ft, r c /r ) are: (a)-(5, 0.40, 
0.86), (b)-(10, 0.74, 0.85) and (c)-(20, 1.56, 0.92). For comparison, the NFW predicted results 
(dotted lines) are also shown with a = a®. The horizontal scales are r/ro and r/r s for the RDH 
and NFW profiles, respectively. 

Fig. 2. — Numerical solutions of gas density and baryon fraction, when the self-gravity of the gas is 
included, for two choices of the central density 5 C = po/p c and the universal baryon fraction fb,BBN- 
The best-fit ft model to each curve is plotted as open squares. The corresponding parameters are 
summarized in Table 1. 

Fig. 3. — An example of the observed and predicted X-ray surface brightness profiles S x for A3158. 
Filled circles: the ROSAT PSPC observed data; Open squares: the best-fit ft model; Dashed line: 
the NFW predicted result; Solid line: the RDH result. Residuals between the best-fit of the RDH 
profile and the data are illustrated in upper panel. 

Fig. 4. — Correlations between the slope and scale parameters in the RDH profile and ft model 
determined from the 45 MME clusters. The solid lines are the best-fit relations. 

Fig. 5. — The same as Fig.4 but for NFW profile and ft model. The dotted line shows the best-fit 
relation by excluding A119, A1367, A1656 and A2255. 

Fig. 6. — Correlations between the density and scale parameters in the RDH and NFW profiles 
derived from the 45 MME clusters. The solid lines represent the best-fit relations. 

Fig. 7. — Central density of dark halos on scales from galaxies to clusters of galaxies. We take 
the data for dwarf galaxies (filled squares), LSB galaxies (open squares) and three distant clusters 
(open circles) from the work of Firmani et al. (2000) directly. Our derived central densities po 
in the RDH profile for 45 MME clusters are shown by filled circles, for which the horizontal axis 
represents velocity dispersion rather than maximum rotation velocity (V m ) as for galaxies. 
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Table 2. Cluster Sample 



cluster 


T (kcV) 


r c (Mpc) 


P 


a 


r s (Mpc) 


Ps (10 4 p c ) 


"0 


A85 


6 74+0.50 

D - ' 4 -0.50 


074+ - 003 

u - u '*-0 .003 


n 5fi2+o o 05 


8.5021°;™ 


S59+0-021 
u.ooa_ Q 021 


1.942lg;|» 


« 795+0- 128 
0. iao_ 12g 


A119 


fi (IT 1 " - 55 


514+0-029 


714+0-025 
u -' 14 -0.025 


14 Q77+1-561 
'-1.561 


O 9 OO+0.480 

o.zoo_ 4g0 


0.038lS;° g 


9 n5fi+°- 288 

a.uoo_ 2gg 


A262 


2 29+ 12 
z -^ a -0.O9 


021+ 002 
u.uzi_ 002 


465+0- 007 
u.^do_ 007 


6.885to;j° 7 


099+ 008 




7 164+ 0121 

/.1D*_ 121 


A401 


10.68l ;i 4 


n 9c -o+0.007 
U " SJJ -0.007 


613+ 005 

U.D±O_ 005 


S 604+°- 100 
8.bU4_ Q 1Q0 


n «^4+0- 030 
U.8t)4_ 030 


n ,-.,+0.102 
u.ooi_ og3 


- 97 +0.064 
8.4J7_ 064 


A426 


7 7, +0.29 
' ■ ' -0.37 


041+ 003 
u.uii_ 0Q3 


n t-o 7 +o.oo7 


7 959+ 105 


192+ 013 
u.iaz_ Q 013 




8 299+° 114 

°-^yy_o.n4 


A478 


7 42+ - 71 

' - 4Z -0.54 


g nl +0.005 
u - 1±i -0.005 


663+ 008 


9.5291°;^ 


406+ 018 
u -™-0 .018 


1.874+g;|g 


9-909t^ 24 4 


A496 




n 042+ 002 


o 539+000 5 


8.152tg;liS 


204+ 014 
u.zu^_ 014 


3.860lg;g 9g 




A644 


7 4 7 +0.32 


n 2n3+ 011 


701+ - 015 
—0.015 


10.153i°; 3 ll 


n 7-JQ+0 055 
U.C3S_ 055 


n 6n9+ 0136 

u.oua_ llg 


Q O55 + 0.070 

a.ooj_ 070 


A754 


I2.85tj; 35 


n 4S1+ - 027 

U.481_ Q 027 


746+ 022 
<™-0.022 


WS89±° f 2 { 


1 646+ 067 

l.DTO_ 067 


n 999 +0.053 
u.^^^_ 046 


9 7Q4+0-275 
a. 1 at_ 275 


A780 


4 49+0- 41 


0.0721°;°°* 


0.6361°;^ 


9 203+ 093 


270+ 012 
u - z ' u -0.012 


9 . 7 .+0.468 
z -^'°-0.446 


9 640+° 115 


A1060 


o 97 +0.11 
° 1 -0.09 


0.0891°;™ 


0.608±g;gli 


8.970t° : i 98 


n oo 7 +0.018 
U.OOI _ 01g 


, , 97 +0.185 
l-l^l -0.178 


8.263±°;i|l 


A1367 


3 99+0- 48 


n 449+0 022 
u.<i<i»_ 022 


n 71 7 +0.024 
u -' ±( -0.024 


19-753tl;| 4 i 


3 969+ 719 
o.aua_ Q 71g 


n22 +0 ' 014 


9 104+0- 292 
a.iut_o 292 


A1651 


7 ic+0.84 


160+ 009 

U.1DU_ oog 


637+ 013 


9 055+0-203 
a.uoo_ 203 


556+ 035 
u.ooo_ 035 


914+ 244 

U.31t_ 21 g 


S fi4q+ 066 
°' uw -0 .066 


A1656 


10.03lg;~ 


o.36ilS;° ° 07 


n fifi^+O-OO 8 
u.ooo_ 00g 


15. 386^1.048 


2 569+0- 235 
i.oua_ 235 


102+0- 035 
uluz -0 .034 


o n/1 +0.100 


A1689 


9.48+ ;f 2 


0.187+^17 


n 7t -n+0.026 
U - ' OZ -0.026 


10 929+0- 456 
iu.y^3_ 456 


0.6801^073 


q 979+0.392 


in 099+ 102 
iu.uaa_ Q 102 


A1795 


7 26+ 51 

'■ ZO -0.40 


n n 9 +0.005 
uilz -0.005 


U.b87_ Q07 


in i4n+ 0133 




i-458i°; 2 S 


10 4m+ 0143 


A2029 


8.22±0-|8 


u - lul -0.006 


633+ 010 


9.12410- 37 


379+0 022 
u -°' -0.022 


oss+0.467 
2.2«8_ 3gl 


9.438l°;i 47 


A2052 


3 30+ 16 
°- ou -0.13 


n n52+°-o 04 

u.uoz_ Q 0Q4 


u.t>»8_ 010 


8.«3tg;l« 


195+ 013 
u.iao_ Q 013 


3 204+ - 633 
°- zu ^-0.603 


8.624l^ 4 7 


A2063 


3 90+ 51 
°- au -0.38 


097+0 06 
uua ' -0.006 


n .-,,,+0.010 

u -°' z -o.oio 


8.073±g;i« 


r, qq 7 +0.020 

u.ooi _ .020 


1.208l°; 32 i 


s.neij;^ 22 


A2142 




u - 13 -0.009 


643+0 09 
u ' u ^°-0.009 


9.200l°;^ 4 4 


699+ 031 
u.uaa_ Q 031 


0.90ll°; 3 ° 2 


9 264+ 048 
a - zu *-0 .048 


A2199 


4 70+° - 13 


067+ 002 
UUD ' -0.002 


6Q0+ 003 
u.ouu_ 003 


« 979+O.O8I 
°- ala _0.081 


307+ - 012 
u.ou(_ 012 




Q OC-q+0.096 

a.ojo_ 096 


A2204 


o ,o + l.(l8 
8 - i8 -1.08 


060+ 006 
u.uou_ 006 


n .07+0.009 
u -°°' -0.009 


cocr + O.lll 

o.ozo_ xll 


226+ 016 


5.985^J 732 


8.962±g;i» 


A2244 


8 47+0- 43 


nl +o.oii 


n 507+O.OI6 


iiq+0.229 
o.ll»_ 22g 


354+ 037 
u.ooi_ 037 


9 qqs+0.690 
^.oao_ 6g7 




A2255 




603+ 043 

U.DUO_ 04 g 


n 01 7 +0.039 
U - 8i ' -0.039 


19 756+ 4 ' 078 
1a. (oo_ 4 07g 


4.57111; 238 


032+ 028 
u.uoz_ 02g 


io.28oig : «i 


A2256 




446+ 018 

u '^™-0.018 


Q 7-Q+O.OI6 
u - ' ' -0.016 


11 941 +°- 546 
n.a<ii_ 54g 


l.768±g;i« 


145+0- 050 
u.liO_ 050 


10.047±°;^ 
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Table 2 — Continued 



cluster T (keV) r c (Mpc) /3 a r s (Mpc) p s (10 4 p c ) a 



A2319 


13.60_ 2 ; 22 


0. 


9 -, +0.014 
ZJ -0.014 


0. 


coo+0.012 
' ooo -0.012 


8.567l°j 2 J 





Q o9+0.063 
aoz -0.063 


0.586tQ^gg 


7 7.9+O.O85 
'• ly:Z -0.085 


A2597 


Q qi +0.27 
J - a -0.22 





067 +0.004 

UD ' -0.004 





660 +o.on 

DDU -0.011 


q 94Q+0.164 
a.zio_ 164 





9, -+0.017 
Z± ' J -0.017 


3 402+ ' 822 


9 640+°- 159 
»-o^u_ 15g 


A3112 


4 69+ ' 27 

u -0.26 





Q54 +0.004 
UJ -0.004 


0. 


t-79+0.007 
O,z -0.007 


a oqn+0.081 

O.OOU g ogl 





9nt -+o.oio 

ZUO -0.010 


4 050+ ' 680 
i.uou_ 672 


« OQQ + 0.100 


A3158 


o oq+1.43 
°- oa -0.98 





245 +0.017 
Z4O -0.017 





64 o+0.019 
uy:O -0.019 


9 485+ ' 521 
».ioo_ 521 





904+°- 108 

au -0.108 


422+ 197 
u.+zz_ Q 172 


a qrq+0.133 
o.ooo_ 133 


A3266 


9 69+ 97 

a - u -0.92 


0. 


c-iQ+0.024 
OiU -0.024 





sn9 +0.021 
8U2 -0.021 


1901 c+0.328 

iz.jio_ 32g 


2 


m 1+0.096 
u±1 -0.096 


129+ 029 


10 311+ ' 234 

1U.01J._q 234 


A3391 


6 90+ 1 ' 47 


0. 


1C7+0.015 
101 -0.015 


0. 


9+O.OI6 
J1 -0.016 


7 991 +0.308 
'■ zzi -0.308 





C79+O.O66 
J,z -0.066 


U.DDO_ 2g4 


fi 709+0.145 

D - ,,;>z -0.145 


A3526 


4 04+° 11 





ni 9+0.001 
ul -0.001 





46 9+0.005 
™ z -0.005 


6 979+ 081 
O-»l»_0 .081 





066+0004 
UDD -0.004 


28.350lt 9 « 7 


7 077+O.IO2 
'-0.102 


A3532 


4 40+ 4 ' 70 


0. 


90,1+0.020 

zo -0.020 





j-qn+0.024 
Jaa -0.024 


9 467+ 742 
».iui -0.742 





99 9+0.148 
yy -0.148 


o.i85l°;f 2 : 


7 0Q7+O.I58 
'• Ja ' -0.158 


A3558 


« Kn+0-50 





lfiCi +0.007 
1DO -0.007 


0. 


54 o+0.007 
o4o -0.007 


7 fil 7+0-098 
' - D± ' -0.098 





560+ 023 
ODU -0.023 


0.699t°;i 2 i 


7 049+O.O54 
'••3^-0.054 


A3562 


6 96+ 1 ' 77 
D - aD -0.95 


0. 


1f1 o+0.006 
±uo -0.006 





48 9+0.006 
^ oz -0.006 


6 7fi q+0.090 
D. <OO_ ogo 





oko+0.020 
• 3o °_0.020 


1 609+ ' 614 

l.DU»_ 424 


6 924+ 062 
u - a ^^_0 .062 


A3571 


o i 9 +0.42 





1 oo+0. 009 
±o ' :> -0.009 


0. 


641+o.on 

D41 -0.011 


9 169+ - 189 

».10»_q lgg 





644 +0.037 
°* i -0.037 


78^+0-148 
U.C80_ 145 


5.00^-0.051 


A3667 


o i , +0.82 
°- il -0.73 





o n q+0.017 
ouo -0.017 





troo+0.012 

°°°-0.012 


8 fiSI +0 269 
8-D81_ 26g 


1 


1 J-9+O.O88 
1JZ -0.088 


n 909+0. 066 
u.zoz_ 063 


7 913+° 147 


A4038 


o oq+1.60 





Q59 +0.004 
uu -0.004 





cc, +0.010 

OO1 -0.010 


7 7/1o +0.134 
'•'^-0.134 





9rn +0.013 
zul -0.013 


9 770+1.747 


'• 80Z -0.115 


A4059 


4.05_°;£ 





ogo +0.005 
ua -0.005 





589 +o.on 

oo3 -0.011 


q-lq+0.156 

o.oto_ 156 





,1 j+0.019 
-0.019 


i.ouu_ 2gl 


01 1 +0.135 
S - Jil -0.135 


AWM7 


3-96±g:iS 





nqq +0.006 
ua -0.006 





607 +0.012 
uul -0.012 


8.863lg;« 4 





oo n +0.023 
■ OOU -0.023 


1 064+ 197 


8 643+°- 147 
o.uio_ 147 


Cygnus 


9 49+ 23 
a -^ a -0.23 





n91 +0.003 
uz± -0.003 





4fis +0.005 
4bS -0.005 


7 17-+O.O66 

'•-"-'^-o.oee 





1n7 +0.007 
lu '-0.007 


26.173l 4 ;^i 


7 fiq7+0-083 
' ■"»' -0.083 


MKW3S 


o +0.16 

^'•"■-o.is 





061+0-002 
uu -0.002 





tro 7 +0.005 
O5 '-0.005 


8-43i±°;iH 





991 +0.012 
zzl -0.012 


9 o m +0.474 

2 -° ui _0.496 


8 521 +°- 123 
o.ozi_ 123 


OphiA 


12.74_°;2° 





17c -+0.006 
11 J -0.006 





644 +0.008 
"H-0.008 


q Ofl 1 +0+38 
'^"-0.138 





614 +0.025 
u± -0.025 


1 Qr:q+0.163 


8 Koq+0-044 
o-oo»_Q 044 


PKS0745 


8-08lg:S 





074 +0.004 
ul -0.004 





614 +0.007 
ul -0.007 


a Q9Q+0.094 
o.ozo_ 094 





9 fi9 +0.013 
zu -0.013 


4 539+ 814 
*' JO3 -0 .769 


9 979+ 0108 


TriaAust 


12.48_ 3 ; 88 





9OK+0.012 
ZOO -0.012 





640 +o.oio 

D ^ u -0.010 


9 436+ - 197 


1 


06g +0.052 
UDa -0.052 


451+ 193 
u -* ol -0 .193 


8.408t°;« 87 



Average 




r (Mpc) 




r c (Mpc) 




r c (Mpc) 




r s (Mpc) 



